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ABSTRACT 

We present MERLIN observations of the continuum (both 1.6 and 5 GHz) and OH maser 
emission towards Arp 220. The correct spatial configuration of the various components of 
the galaxy is revealed. In the eastern component the masers are shown to be generally coin- 
cident with the larger-scale continuum emission; in the west, the masers and continuum do 
not generally arise from the same location. A velocity gradient ((0.32±0.03)kms^^pc^^) is 
found in the eastern nuclear region on MERLIN scales; this gradient is three times smaller 
than that seen in HI and implies that the OH gas lies i nside the HI. A re-analysis of previ- 
ously presented global VLBI data ( Lonsd ale et alJl998l) reveals a very high velocity gradient 
((18.67±0.12)kms^^pc^^) in one component, possibly the site of a heavily obscured AGN. 

Key words: masers - galaxies:individual:Arp220 - galaxies:active - radio contin- 
uum:galaxies - radio lines:galaxies. 



1 INTRODUCTION 

Ultra-luminous infrared galaxies (ULIRGs) are so named because 
their infrared luminosity is greater than the luminosity in all the 
' other wavelengths combined (L/_R > lO^^L©). The source of such 
copious emission is thought to be the enormous amount of molec- 
ular dust in their central regions, which obscures visible and other 
wavelengths and produces infrared light through heating. The dust 
grains are excellent absorbers of short wavelength radiation emit- 
ted from newly formed OB stars, so a high infrared luminosity is 
a tracer of star formation. In many cases however active galac- 
tic nuclei and Seyfert-like nuclei can play an important role. A 
review of the properties of this kind of galaxy can be found in 
ISanders & MirabelHl99d) . 

One of the most interesting properties of ULIRGs is that they 
tend to be merging systems, as it takes a strong merger of dust 
rich galaxies to generate the observed high density star formation. 
The proportion of the mergers among ULIRGs gets larger when the 
infrared luminosity gets higher, and is 100% in luminosities above 
10^^1/0. (See Sanders & Mirabel ( 1996) and references therein). 
Many of them also have ve ry powerful OH maser emission (e.g. 
iDarling & Giovanellil i20o3) ) in which the luminosity exceeds that 
of known galactic masers by a factor of at least 10^, hence the term 
"megamaser". 

Arp 220 is one of the most luminous infrared galaxies with 
LiR = 1.4 X 10^^ L0 iSoiferetal]ll98"i) and provides a superb 
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example of the starburst phenomenon, demo nstrated graphicall y 
by the discovery of several radio supernovae ^Smith et al."l99d). 
Starburst p henomena are used b y [Shiova, Trentham & Tanigushj 
(2001) and llwasawa et al] ^2001) to explain the X-ray properties 
of Arp 220, although the existence of a hea vily obscured AGN is 
not ruled out. Recent CHANDRA results iClements et alJl2002h 
show that the hard X-Ray emission in Arp 220 is confined in a 
sub-kiloparsec scale region, in contrast to other starburst galaxies 
and its spectrum shows that the hard X-rays are more likely to be 
produced by one or more low luminosity, heavily obscured AGN, 
X-Ray binaries or Ultra Luminous X-Ray sources rather than the 
supernovae. Therefore the co-existence of such structures is still 
plausible. 

The morphology of Arp 220 is also an open issue. At op- 
tical I'Scoville e t alj|l993) as well as at radio wavelengths (eg. 
P^orris 1988) it exhibits a double structure with tidal tails and dust 
lanes which leads to the assumption that it is the result of a re- 
cent merging phenomenon. Studies of CO, HI and the OH maser 
emission enable the determination of its velocity structures. On 
balance, the current data appear to favor a model of two coun- 
terrotating disks as parts o f a larger scale disk (Sakamoto e t all 
[l999; Mu ndell. Permit & P edlar 200lJ), but a single warped cir- 
cumnuclear disk model can also be fit to the data lEckart & Downe j 
2001). In the case of the related galaxy IIIZw35 the existence of ev- 
idence for such a disk, where the OH maser is located, was recently 
revealed (j'hilstrom et al. 2001) so a search for a similar structure 
in Arp 220 would be useful. 

The very powerful maser emission is perhaps the most 
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striking property of Arp 220. It was first discovered by 
iBaan. Wo od & Haschick 1 198^ while searching for OH absorption 
in galaxies with known HI absorption. The masers were th ought 
to ari se in regions extending to several hundreds of parsecs jBaanI 
I l98^. However, hig her resolution studies iDiam oiid et alJll98^ 
iLonsdale et aljl99a) revealed a dual component structure: there is 
a diffuse component which consists of approximately one third of 
the 1667 MHz OH emission plus almost all of the 1665 MHz OH 
emission with a physical extent of some hundred parsecs as sug- 
gested by the Baan model; but there is also a compact component 
consisting only of the 1667 MHz line with a physical extent of just 
a few parsecs. These compact masers show high velo city widths, 
even in this very small physical extent iLonsdale et 31119981) . This 
means that either they come from a region with a volume of a 
few pc'^ or they represent long, thin filaments with their orienta- 
tion along the line of sight. In the first case there has to be a hidden 
infrared source in that region to provide the essential pumping pho- 
tons, or the optical depth of this region has to be large enough to 
trap the infrared photons, if the pumping is indeed radiative and 
not coUisional due to shock fronts from the nearby supemovae. 
If filamental structure was present to any great extent we would 
have to explain why the filaments' orientation is onl y along the line 
of sig ht, so that we cannot see elongated structures iLonsdale et all 
I2OOII) . Probably we are only seing filaments which lie on the line 
of sight due to the gain path needed for a maser to be detected. 

In this paper we present new MERLIN OH spectral line and 
continuum data, and compare it with previously p resented global 
VLBI data jSmith et all 1991 Lonsdale et alll998h . 

2 OBSERVATIONS AND DATA REDUCTION 

In this section we describe the MERLIN and global VLBI observa- 
tions of both line and continuum emission 

2.1 MERLIN Observations 

2.1.1 1.6 GHz Spectral Line Data 

Data were taken on 23 April 1998 with seven MERLIN antennas 
including the 76 m Lovell telescope observing with left-hand cir- 
cular polarization only. The dataset included observations of the 
target source Arp 220 (spending 9 hours and 20 minutes on source) 
and the phase reference 1511-1-238, as well as two scans on the 
point source calibrator 2134-1-004 and one scan on the flux calibra- 
tor 3C286. The observations covered a bandwidth of 4 MHz and 
were correlated to produce 256 channels, each with a bandwidth of 
15.625 KHz. 

The initial editing and gain elevation corrections, as well as 
preliminary amplitude calibration were made using MERLIN spe- 
cific software. Then, the data were transferred into NRAO's Astro- 
nomical Image Processing System (AIPS) for bandpass calibration 
and phase corrections. Phase referencing was performed using the 
data on 15 1 1-1-238. The data were then self-calibrated starting from 
a point source model and the solutions were applied to the Arp 220 
data. After that line free channels were averaged and an image of 
this continuum was produced in 'difmap'. There may be modest 
contamination of the continuum image by weak 1665 MHz OH 
emission but tests have shown this to be minimal. We attempted to 
image the 1665 MHz emission but it was too diffuse even on MER- 
LIN scales; it was detected only on short baselines, shorter than 
0.35 A/A. Several iterations of self-calibration were performed on 



the continuum data and the solutions obtained were then applied to 
the multi-frequency dataset and the data averaged for 60 seconds. 
The final spectral line dataset was generated by subtracting the con- 
tinuum emission using UVLIN and then averaging every two chan- 
nels together, resulting in a velocity resolution of 5.6km s^^. 

2.1.2 1.6 GHz Continuum Data 

Continuum data were taken on 23 June 1998 with six MERLIN an- 
tennas using right hand circular polarization. The dataset includes 
observations of the target source Arp 220 (spending 4 hours on 
source) and the phase reference 151 1-1-238, as well as two scans on 
the point source calibrator 2134-1-004 and two scans on the flux cal- 
ibrator 3C286. 16 channels were generated by the correlator, each 
with a bandwidth of 1 MHz, the central channel (8) being set at a 
frequency of 1666.00 MHz. 

Local MERLIN software was used to perform the gain eleva- 
tion corrections and preliminary amplitude calibration and editing. 
Bandpass calibration was performed and AIPS was used for fur- 
ther calibration. After flux calibration, the phase reference source 
was used for phase calibration, and three rounds of self-calibration 
were performed on 1511-1-238 starting from a point-source model. 
The solutions were copied and applied to all the data, channels 2 to 
15 were averaged and a first image of Arp 220 was produced. Start- 
ing from this image, we performed four iterations of self calibration 
on the Arp 220 data before creating the final image. 

2.1.3 5 GHz Continuum Data 

5 GHz data were taken using 6 MERLIN antennas in full polar- 
ization on 29 May 2000. The dataset consisted of observations 
of Arp 220 (spending 7 hours on source) and the phase reference 
1551-1-239 as well as five scans on the point source 2134-1-004 and 
two scans on the flux calibrator 3C286. Again, 16 channels were 
produced by the correlator each with a bandwidth of 1 MHz, the 
central channel being set at a frequency of 4.995 GHz. 

Standard procedures were followed for the data reduction, the 
phase corrections from the phase reference source were applied to 
the Arp 220 data and after self-calibration the final image was pro- 
duced using Stokes I. 

2.2 VLBI Observations 

The VLBI data used for comparison ar e those discussed in ou r pre- 
vious papers, Smith etal. (1998) and lLonsdale et all il998h . The 
data were taken in November 1994 and discussion of the calibra- 
tion is contained with the papers. 

3 RESULTS 
3.1 Continuum 

3.1.1 MERLIN 

The 1.6 GHz (June 1998) MERLIN continuum image of 
Arp 220 is shown in Figure 1: A double component structure 
is present, with the western and the eastern components sepa- 
rated by ~1 arcsec (~370pc if we assume a distance of 76Mpc; 
H()=75 kms~^ Mpc~^) and at the same positions as the 1.3 mm 
( Sakamoto et al. 1999) and the 4.83 GHz (Baan & Haschic Jl99l 
this paper) continuum peaks. The integrated flux densities are 
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Figure 1. The 1.6 GHz image of Arp 220 from the MERLIN continuum 
data. The peak flux is 44.604 mJy beam" ^ and the contour levels are at 1, 
2, 4, 8, 16 and 32mjy. The beam size is shown at the bottom left corner 
of the image and is 209.97 X 167.48 mas (77.4x61.7 pc) and the off source 
noise level is 055 mjy beam~^ . 



(91.9±0.9) mJy and (1 1 1.2±0.8) mJy for the eastern and the west- 
em components respectively. We should note here that the com- 
ponent boundaries are somewhat arbitrary, so the overall diffuse 
emission plays a role in deriving these numbers: the values given 
are lower limits. The total integrated flux density of the galaxy at 
1.6 GHz is (280±40)mJy, consistent with that derived at 1.4 GHz 
bv lMundelletalJl200lh (285.4±14.3 mJy). 

In Figure 1 we can see that the eastern component is elon- 
gated in a northeast - southwest direction and a 'spur' of emis- 
sion is visible on the northwest of the western component. This 
is also detected by Mundell et al. (20()j) at 1.4 GHz and by 
iBaan & Haschickl fl995l) in the 4.8 GHz continuum. A compo- 
nent southeast of the western component is also detected as part 
of a more extended emission p attern, possibly part of th e extended 
emission detected by the VLA (Baan & Haschick 1984). The peak 
flux of this region is 8.72mJybeam^^, whereas the off-source 
noise level is 0.55 mJy beam^^. 

In Figure 2 we can see the 5 GHz image of Arp 220 where 
again the double component structure is present, although more ex- 
tended structure is not detected. The peaks of the two components 
are separated by (1.066±0.004) arcsec and their integrated fluxes 
are (33.6±1.0) mJy and (60.7±0.8) mJy for the east and west re- 
spectively, while the off-source noise level is 0.4 mJy beam" ^. If 
we compare these fluxes with these of the 1.6 GHz data we can de- 
rive the spectral indices of these two continuum components to be: 

= -0.90 ±0.03 
Q„ = -0.542 ±0.013 

The errors shown here are statistical. There might be systematic ef- 
fects that cause larger errors. The spectral indices derived here are 
similar to the spectral index (for the whole nucleus) assumed by 
iNorrisliTgsl (-0.6 to -0.9). This is a characteristic of synchrotron 
radiation (see next paragraph) though it cannot give us a safe indi- 
cation of the exact origin of the emission or its power source. 
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Figure 2. The 5 GHz image of Arp 220 from the MERLIN continuum data. 
The peak flux is 34.6 mJy beam" ^ and the contour levels are at 1, 2, 4, 
8, 16, and 32mJy, as in Figure 1. The beam size is shown at the bottom 
left comer of the image and is 178.42x95.16 mas (65.7x35.1 pc). The off- 
source noise level is 0.4 mJy beam^ ^ . 



3.1.2 Global VLBI 

The results fr om the continuum V LBI dataset have been previously 
presented by ISmith et'al] il998h . where the discovery of proba- 
ble radio supemovae (RSNe) is reported. There are about a dozen 
RSNe (November 1994) in the western nuclear region and possibly 
one or two in the eastern, without any other obvious compact con- 
tinuum component. (Figures la and lb respectively in Smith et al] 
( 1998)) The sum of the flux densities of the observed supernovae 
in the western region is 7.87 mJy and in the eastern is 0.78 mJy in- 
cluding the marginally detected ones. If we compare these with the 
flux densities observed with MERLIN it is safe to assume that the 
MERLIN and VLBI continuum detections are the results of differ- 
ent processes; the large-scale MERLIN continuum emission proba- 
bly comes from synchrotron acceleration of free electrons (possibly 
generated in the supernova explosions), while the VLBI emission 
probably arises from the interaction between the explosion material 
and the stellar wind generated before the star exploded. 



3.2 Spectral Line 

3.2.1 MERLIN 

The naturally weighted velocity maps of Arp 220 (after continuum 
subtraction) can be seen in Figure 3. There are three components 
present, one in the east and two in the west (northwest and south- 
west). The eastern component is more redshifted than those in the 
west, suggesting an overall velocity gradient along the east-west 
axis. This is i n agreement with the suggested presence of an overall 
rotating disk ('Scovill e. Yun & Brvanal997tlSakamoto et alll999t 
Mundell et al. 200 1 ). The off-source noise level is 3.2 mJy beam" ^ . 

In Figure 4 we can see the spectral line emission (naturally 
weighted and averaged over the velocity range of the OH emission) 
superimposed on the continuum emission. This continuum image 
(in greyscale) is not the same as that in Figure 1 , but was generated 
from the continuum portion of the spectral line dataset. The error in 
relative position between the continuum and the spectral line image 
is less than a milliarcsecond but the noise in this continuum image 
is significantly greater than that obtained from the wideband con- 
tinuum data (0.83 mJy beam"^). The eastern portion of the spectral 
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Figure 3. The natural weighted images of Arp 220 from the MERLIN spectral line data. The peak integrated flux is 166.73 mjy beam ^ . The beam size is 
shown at the bottom right comer of the first image and is 229.96 X 182.75 mas (84.7x67.3 pc). The noise level of each image is 3.2 mJy beam~^. 



line emission is seen to be at the same position as the continuum, 
but the peak of the western component of the continuum emission 
lies between the two regions of spectral line emission. The north- 
western line component is at the same position as the compact 
VLB! component and the southwestern emission corresponds to 
the previously observed diffuse component (Lonsdale et al. 1998). 
The eastern component spectrum is broad and has a double peak, 
which implies that regions with different velocities are superim- 
posed (visible with the VLBI) to produce the component detected 
with MERLIN; their spectra are blended together. 

3.2.2 Global VLBI 

iLonsdale et alj il998h presented the results of the VLBI spectral 
line dataset, where they find both a diffuse and a compact com- 
ponent in each of the eastern and the western regions. In Figure 5 
we can see the spectral line VLBI image of Arp 220 (in contours) 
with respect to the continuum emission (in greyscale). Both images 
are convolved with a taper of 10 MA to include both regions in one 
map, so they are of degraded spatial resolution. In the western re- 
gion we can detect the maser in regions where we have continuum 
emission due to the supernovae, but the majority of the compact 
maser emission is detected in regions north and south of the super- 
novae, where no compact continuum is detected. In the case of the 
eastern region there is no clear compact continuum emission, only 
hints of possible radio supernovae; the maser does not follow the 



pattern seen in the west. Close-ups of the two nuc lear regions are 
shown in contours (the same as in figures 1 and 2 of lLonsdale et alj 
with full resolution. (Beam size 9.8x2.5 mas). Spectra of 
different regions are also shown. One interesting point to note is 
the fourth spectrum (lower left in Figure 5), which shows absorp- 
tion. We should note that there are problems with imaging diffuse 
emission with the global VLBI due to the lack of short uv spac- 
ings, we are however confident of some absorption being present. 
This suggests that the OH radical is present even in regions with 
no maser emission. We also detect part of the 1665 MHz emission, 
clearly seen in the two diffuse regions (spectra 5 and 8) and in the 
eastern compact region (spectra 1 and 2). To determine the origin 
of this emission we used tighter boxes to determine the spectra; 
the 1665 MHz emission then fades dramatically with respect to the 
1667 MHz emission. This suggests that it is diffuse and comes from 
the overall area, a result first reported in iLonsdale et aL ( .1998f) . 

4 DISCUSSION 

4.1 Positions of the Masers 

Any discussion of positions must be prefaced by the comment that 
since we performed self calibration, information about the absolute 
positions of the components is lost. The conclusions we make be- 
low are based on the relative positions of the components. We have 
extremely accurate relative positions (sub-mas) since we are able 
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(5RA (arcsec) 


(5DEC (arcsec) 


NW maser W cont 


MERLIN 


0.011 ±0.006 


0.191 ±0.004 




VLBI 


0.0024±0.0007 


0.1394±0.0007 


NW maser SW maser 


MERLIN 


-(0.012±0.011) 


0.311 ±0.007 




VLBI 


-(0.0335±0.0007) 


0.3032±0.0016 


W cont SW maser 


MERLIN 


-(0.022±0.008) 


0.12±0.02 




VLBI 


-(0.0359±0.0012) 


0.1638±0.0017 


NW maser E maser 


MERLIN 


-(0.993±0.006) 


0.339±0.004 




VLBI 


-(1.0710±0.0002) 


0.3252±0.0006 



Table 1. The positional differences between the individual components in both the MERLIN and VLBI spectral line and continuum emission. 



to obtain low sensitivity images of the continuum emission from 
the spectral line datasets; these can be used to tie down the position 
of the high sensitivity continuum image with respect to the maser 
positions. 

It is clear that in both MERLIN and VLBI maps the contin- 
uum emission does not coincide with the maser peaks. Table 1 lists 
the measured offsets between various maser and continuum com- 
ponents in the MERLIN and VLBI images. (We chose the brightest 
supernova as the location of the VLBI western continuum peak and 
the brightest maser component (area of spectrum 2) as the peak of 
the spectral line emission in the east). 

From this table we can assume that the northwestern MER- 
LIN component of the maser emission is at the same position as the 
compact northwestern VLBI maser component and the southwest- 
em maser MERLIN component is at the same position as the dif- 
fuse southwestern VLBI maser component. In this case, the bright- 
est radio supernova is on the northern part of the western MER- 
LIN continuum emission, ~0.05 arcsec offset from its center. In 
the eastern region, the VLBI maser emission is situated on the bor- 
der of the MERLIN maser emission. Since the MERLIN maser and 
bulk of the continuum coincide within ^--^0.07 arcsec, and the possi- 
ble supemovae are situated in the area between the two VLBI maser 
components ( Smith et al. 1998), they lie in the area of the MERLIN 
continuum emission. 

Unfortunately, this registration of the different components 
in Arp 220 r e quires a revision of the alignment suggested by 
iMundell et alj J200ll) (see their figure 8). We are confident that 
our alignment is correct, as we extracted both the continuum and 
line images of Figure 4 from the same dataset performing the same 
self calibration. This may, to some degree, affect their conclusions 
regarding the location of the supernovae within their tilted disc 
model. 



4.1.1 Eastern Region 

If we assume that the MERLIN 1.6 GHz continuum emiss ion coin - 
cides with the 1.4 GHz continuum emission (Mundell et jlj20^^ 
then all 1.6 GHz c ontinuum, 1.4 GHz, 1.3 mm iSakamoto et alj 
Il999h and 4.8 GHz feaan & Haschicklll995h come from the same 
area, as does the HI absorpt ion jMundell et all200lh and CO emis- 
sion iSakamoto eta lll999h . Weak, tentatively detected supemovae 
are also situated there. 

The MERLIN maser emission coincides with the continuum 
in the eastem component. The VLBI maser emission lies in re- 
gions with no obvious compact continuum background; it does not 
cover all the area covered by the MERLIN maser but is instead 



situated on the edge of this area having a face-on disk-like appear- 
ance. It is probably saturated; the lower limit of the amplification 
factor is 1 10 for the northem and 65 for the southem region (for the 
1667 MHz line). However, the OH radical is present throughout the 
whole area, as exemplified by the detection of OH absorption. 

4.1.2 Western Region 

In the westem region, the MERLIN 1.6 GHz, the 1.4 GHz, 4.8 GHz 
and 1.3 mm continuum emission come from the same area as the 
CO emission and HI absorption. The strong supemovae are situated 
there, close to the northern border, probably assisting the young OB 
stars in heating the dust grains. The position of the OH maser emis- 
sion is however different. Both MERLIN and VLBI demonstrate 
that it arises from regions north and south of the continuum. The 
lower limit for the amplification factor (derived from the VLBI 
data) is 125 for the compact northern region, 45 for the southem 
diffuse region and 90 for the compact region on its southwest bor- 
der, so that the two compact regions are very likely to be saturated. 
OH is present at some level throughout the westem region; we can 
see some maser emission lying in front of two supernovae northeast 
and northwest from the diffuse maser region. However, the ampli- 
fication factors at the line peaks in front of the supemovae are very 
small, 4.6±0.7 and 3.3±0.3 respectively. This means that they lie 
in the area of diffuse, low-amplifying OH gas, which we only can 
detect in those regions with significant continuum radiation. Obser- 
vations with better sensitivity will be able to detect more sources of 
this kind and so study the properties of the diffuse OH in some 
detail. 



4.2 Velocity Gradiens 

Figure 6 shows the moment map taken from the MERLIN data. 
There is a clear velocity gradient of (0.32±0. 03)kms~^pc~ ^ in the 
eastem region, at a position angle of 35.5° . iMundell et all i200lh 
find (1.01±0.02)kms"^pc"^ for the HI at a PA of ~55°. If we 
assume that the OH arises from a circular rotating disc, its radius 
is (79±5)pc and its dynamical mass is (12±3)x IO^Mq. This is 
smaller by a factor of 90 than the mass of 1.1 x IO^Mq calculated 
from HI absorption. This large difference may be explained if the 
OH gas lies in the inner part of the disk and the HI and CO arise in 
the outer part. This disc may be warped, explaining the difference 
in p osition angle between the OH and HI gradients. 

IMundell et all <200ll) and lSakamoto et alj <1999l) find a clear 
velocity gradient in HI and CO in the westem region of Arp 220, 
in a roughly horizontal direction. In the MERLIN image of the OH 
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Figure 6. The moment map of Arp 220 from the MERLIN spectral Hne 
data. The intensity represents the flux density and the colors the velocity 
field, the reddest being the most redshifted. The contours represent the ve- 
locity field and they range from 5305 to 5475 km s~^ with spacings of 10 
kms~^. The area selected to plot the velocity field is where the integrated 
flux is more than 15 mJy beam^^. 
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Figure 7. The moment maps of the western region of Arp 220 from 
the VLBI dataset. The left picture (a) represents the northern compact 
region and the right (b) the southern diffuse region. The contours for 
both pictures are 12.5kms~^ apart, representing a velocity field from 
5250 to 5612.5 kms~^. Only data whose integrated flux is greater than 
1.5mjybeam~^ are selected. 



(Figure 6) we can see some velocity structure in the southwest- 
em component, although a gradient is not clearly visible and in the 
northwestern component no gradient is detectable. However, the 
VLBI moment maps of Figures 7a and 7b (of the northeiTi and 
southern region respectively) show that the southern component 
has a roughly E-W velocity gradient of (18.67±0. 12)kms~^pc~^ 
(Figure 8). The region from which this gradient arises is on the 
southern border of the region that contains the HI and CO gradi- 
ent and its diameter is (5.9±0.5) pc. If this gas was in rotation, the 
enclosed mass would be (1.7±0.4)x 1O^M0. The velocity gradient 
we measure is ~20 times greater and the mass is ~ 10 times smaller 
than that measured using HI absorption (0.83±0.02)kms~^pc~^ 




Figure 8. The velocity gradient of the southwestern region of Arp 220. 
The points from which the velocities were taken lie on the yellow line of 
Figure 7b with a position angle of 85.2°and no smoothing has been ap- 
plied. The velocity gradient was derived using the least square method and 
is(18.67±0.12)kms-lpc-i. 

and 1.7xlO*M0 respectively. A plausible explanation is that the 
OH is confined in a region of much smaller physical extent than 
the HI, in the center of the molecular disk. However, in comparison 
to the eastern region the mass is concentrated in the centre of the 
disk and that gives rise to the very high velocity gradient. Such a 
disk is however not well confined, so we cannot derive reliable re- 
sults about the mass or determine how the density decreases with 
distance. 

The northwestern component shows a very strange velocity 
structure. On the MERLIN scales no obvious structure is detected, 
whereas the VLBI structure can be seen in Figure 9. If we accept 
that the systemic velocity of this region is ~5355kms^^ (ie. the 
velocity of the central part of Figure 7a), then the region just out- 
side and on both sides of the centre is blueshifted and it becomes 
redshifted as we go further out. Two possible explanations are pro- 
posed to explain this structure, the first being that the masing gas 
experiences spiral outflow within a cone, thus exhibiting an expand- 
ing helix structure (Figure 10). The problem with this configuration 
is that the velocity structure is symmetrical, so that one helix would 
be a mirror image of the other. Alternatively, there might be out- 
going and ingoing spherical shock-waves producing the redshifted 
and blueshifted velocities (Figure 11). The observation of an elon- 
gated structure would be due to beaming effects. 

To distinguish between these two proposed or other possible 
structures and investigate the nature of the central region higher 
spatial resolution observations have to be obtained, or a better 
signal-to-noise ratio has to be achieved. We will then be able to 
examine the velocity structure further away from the central point 
to see if it continues to be redshifted or instead becomes blueshifted 
again, something which would favor the outflow scenario. 

4.3 Is There Still Room for an AGN? 

When OH megamasers were first discovered, the question of 
whether a starburst or an AGN were providing t he essential energy 
was raised. The discovery of RSNe in Arp 220 iSmithetalJl998l) 
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Figure 9. The velocity structure of the northwestern region of Arp 220. The 
different velocities are taken from points along the yellow line of Figure 7a 
which lies along the major axis of the emission. No smoothing has been 
applied. 
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Figure 10. The expanding helix model as an explanation for the VLBI 
northwestern velocity structure. Note that the diagram and the sketch of 
the region are not in scale. 



answers this question in favour of the starburst. Luminous radio 
supemovae or nested supernovae remnants have since been tenta- 
tively detected in a number of OHMs (M rk 273 ^arilli & Tavlor 
12000) . I RAS 17028-0014 iMomiian et al. 2003 in Dreparation), 
Mrk 23l lTavloret"aLl < 19991) and IIIZw35 (Philstrom et aL.2001.) ) 
linking the megamaser phenomenon with a nuclear starburst. The 
infrared luminosity is thought to arise from the vast amount of 




Figure 11. The double shock front scenario as an explanation of the north- 
western area VLBI velocity structure. Bold lined indicate the regions which 
the masers we observe come from. The diagram and the sketch of the region 
are not in scale as in Figure 10. 



dust resulting from merging phenomena heated from young mas- 
sive stars, whereas the radio continuum comes from relativistic 
electrons (partially generated from young stars and partially from 
supernova explosions) being accelerated by magnetic fields, so 
there is a correlation between the radio and the infrared fluxes 
<Bressan. Silva & Granatcl2002h . 

An interesting common feature common of OH megamasers 
is that the emission comes from regions with clear velocity gra- 
dients, interpreted as molecula r discs or torii. Thi s is evident in 
Arp 2 20 (this paper), IIIZw35 tPhilstrom et alj200lt) . I RAS 17028- 
0014 jMomiian etai] 2003 in preparationh. Mrk 273 jYates et alJ 



12000) and Mrk 231 iCarilli. Wrobel & Ulvesta(l"l998Y The for- 
mation of such a disc (or torus) is modeled by Barnes 1 2002) for 
merging gas-rich galaxies and a similar tor us seems to be in agree - 
ment with high-resolution data of IIIZw35 <Philstr6m et all200lh . 
Given the high column densities of these regi ons (8x10^'^ T,(K)< 
NiY < 2.4x10^" Ts(K)cm~2 for Arp 220 HI <Mundell et al J200lh ) 
they can obscure low luminosity AGN. Indirect evidence in the 
form of a radio jet i n Mrk 231 suggests that it hosts an AGN 
<UlvestadetalJll999h . Compact hard X-ray sources i n the nu- 
clei o f Arp 220 jClements et al] Eoo3) and Mrk 273 iXia et alJ 
l2002ft suggest that there is an active nucleus in those galaxies. 
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The existence of an AGN in OHM galaxies may well be a ran- 
dom phenomenon not associated with their masing properties, 
however statistical analysis of masing and non-masing LIRGs 
by Darling & Giovanelli (2002) reveals that LIRGs that host OH 
megamasers tend to be of infrared excess compared to those that 
don't. This means that there is another process taking place which 
contributes to the infrared flux and not the radio continuum. This 
can be a hidden AGN or a very recent starburst, which has not 
yet produced supernovae and relativistic electrons, OH megamasers 
can work as tracer for this kind of phenomena. To distinguish be- 
tween the two possibilities one needs to model violent and sudden 
starbursts to check if the time scales until it starts contributing to 
the radio flux can explain the ratio of OHMs among LIRGs. Al- 
ternatively, a large number of OHMs and host galaxies should be 
observed with high sensitivity to determine if they host compact 
regions, candidates for hosting an AGN. 



5 SUMMARY 

Imaging of both the continuum and spectral line (OH) emission of 
Arp 220 using the same dataset has revealed the correct relative 
positions of its components. In the eastern radio nucleus the con- 
tinuum and spectral emission coincide, whereas in the western the 
continuum is situated between the two spectral line emission peaks. 
On VLBI scales the continuum emission is in the form of several 
point sources, probably radio supernovae JSmithetalll998t) in the 
same position as the MERLIN emission. 

On large (MERLIN) scales the maser emission comes from 
three main regions, two western and one eastern. The eastern re- 
gion has a very clear velocity gradient, sho wing rotation of the nu- 
clear disc, if we foll ow the standard model iSakamoto et alll999t 
iMundell et aT]l2 001^. The velocity gradient we measure is sig nifi- 
cantly smaller than that found using HI absorption I Mundell et alj 
I2OOII) . This may mean that the OH rotates inside the HI on large 
scales. 

On smaller (VLBI) scales the OH maser emission has pecu- 
liar characteristics. In both nuclear regions the 1667 MHz line is 
emitted from regions of very small physical extent, whereas the 
1665 MHz line seems to be of a more extended nature, coming 
from the overall region. Even on MERLIN scales (several hundred 
milliarcseconds) its location and structure could not be reliably de- 
termined. The most luminous maser is situated in the northwestern 
region and its velocity structure (Figure 9) implies that it either 
has a cone-like structure, which could be the result of molecular 
outflow from the central region or that it is produced by a pair of 
ingoing and outgoing shock fronts. 

The issue of compact continuum structures in Arp 220 still 
remains open. Although a starburst, as exemplified by the exis- 
tence of radio supernovae, can explain Arp 220s bolometric lu- 
minosity l'Sm ithetalJll99 ^. the hard X-Ray emission spectrum 
and spatial distribution show that they are unlikely to be its main 
power source and other alternatives are proposed, for example sev- 
eral Ultra-L uminous X-Ray sourc es (ULX) or one (or several) low- 
mass AGN jClements et alj20o3) . The very high velocity gradient 
found in the southwestern maser (Figure 7b) marks it as a possible 
candidate to host a ULX or a low mass (~ 1.7 X 10^ A/q) AGN. 
The sites of the radio supernovae are of very small extent and show 
very faint maser emission, so it is not feasible to study velocity gra- 
dients in these regions. Long time monitoring and analysis of their 
light curves is needed for their properties to be revealed in detail. 
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